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Abstract 

We present a detailed study of the entanglement dynamics of a two-qubit system coupled to 
independent non-Markovian environments, employing hierarchy equations. This recently developed 
theoretical treatment can conveniently solve non-Markovian problems and take into consideration 
the correlation between the system and bath in an initial state. We concentrate on calculating 
the death and rebirth time points of the entanglement to obtain a general view of the concurrence 
curve and explore the behavior of entanglement dynamics with respect to the coupling strength, 
the characteristic frequency of the noise bath and the environment temperature. 
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I. INTRODUCTION 



Quantum information, including quantum entanglement and quantum dissonance [l|-|3j, 
has been a popular area of research in recent years. In this field, entanglement undoubtedly 



plays a central 
dense coding 



role as an important resource in quantum computation 



4|, teleportation 



6] and cryptography [7|, |8[ . There have been many meaningful works on the 
dynamics of entanglement in multiple qubits, especially two qubits, interacting with different 
kinds of environments and we now 
as the entanglement sudden death 

However, environments with different properties will have different effects on the entan- 
glement dynamics. Recentl y q uantum systems in a non-Markovian environment have been a 



mow some important features of the entanglement such 
9-12] and birth Q. 



subject of great interest 
this situation, such as Ref 



Theoretical treatments have been developed to deal with 



161 ] and 20]. Of particular interest to us is the hierarchy 



equations approach employed by Dijkstra and Tanimura (l4j. This treatment is deduced 
using the influence functional method of Feynman and Vernon without the limitation of 
perturbative, Markovian or rotating wave approximations. It can easily take into account 
the effect of the system-bath coupling on the dynamics of the entanglement for any initial 
conditions. In their work, they discussed the entanglement evolution of two qubits inter- 
acting with a quantum-mechanical bath and then compared this hierarchy method with the 
Redfield equation. It is found that the result of the full calculation markedly differs from 
the Redfield predictions. In the present work we continue this meaningful work. 

We first introduce the method of hierarchy equations of motion developed by Y. Tanimura 



and coworkers in Ref. 



2lM23l|. Then we employ this method to carefully calculate the en- 



tanglement dynamics of a two-qubit system interacting with a non-Markovian environment, 
determining the influence of the strength of the system-bath interaction, the characteris- 
tic frequency of the bath and the environment temperature on the time evolution of the 
entanglement, especially on the sudden death and sudden birth time points. 

The entanglement of the two qubits should be measured using Wootters' concurrence 



C(p) = max(0, yAi 



A 4 ), 



(1) 



2 



and 

P = PAB(tf®<Ty)p* AB ((7y®<7y), (2) 

where pab is the density matrix of system AB, p* AB denotes the complex conjugation of 
Pab and is the Pauli matrix. \ are the eigenvalues of p (Ai should be the largest 

eigenvalue). 



II. THE MODEL AND THEORY 



For simplicity, we set h = 1 throughout this report. We assume that the two qubits are 
coupled independently to two identical baths with the same strength. The baths have a 
characteristic frequency 7, and large and small 7 indicate fast and slow baths respectively. 
The energy gap of the qubits is e and the qubits are coupled by an interaction (. Therefore, 
we can write the standard system Hamilton as [3] 

H s = e(a\ai + a\a 2 ) + ((a\ + (Xi)(a\ + a 2 ). (3) 

The subscripts 1 and 2 represent the two qubits and a) and a are the creation and annihilation 
operators. We choose a model that adequately represents the environment as a set of 



oscillators with a coupling linear in the oscillator coordinates, having the form 26] 

H * = + ^-M) ( 4 ) 
3 3 

and 

H S b ^ c mj{a ] m + a m )xj, (5) 

m,j 

where Hb is the Hamilton for the bath and H$b is the interaction item of the system and 
bath, and Uj are the coordinate, momentum, mass, and frequency of the jth 

harmonic oscillator, respectively. C m j is the strength of coupling of the mth qubit to the jth 
oscillator. To obtain complete information about the effect of the environment, we introduce 
the spectral density J m (co), which is defined as 

Jnii^) = I ^(C^/maUaWu ~ U a ), (6) 



3 



and in this report we assume that 

J m (u) = uj^/iu 2 + 7 2 ) (7) 



which comes from the Lorentzian cutoff [22j, and is the same for both qubits. The the 
hierarchy equations approach then gives the equation of motion for the system density 
matrix, which has the form Q, 



25j: 



1 /i\ 2 K 2 _ , K 

dp n {t) _ 



£ + ^2 n mklk)pn(t) ~ — - ~ —Wm, [V m , Pn{ 

m=l k=0 m=l k=0 ^ k 

2 K 2 K 

^Pn mk +l(t)) — i y^y^n m k(CkV m Pn mk -l{t) — C* k p nmk -i(t)V m ) , (8) 
m=l k=0 m=l k=0 

where Cp = [Hs,p], 7o = 7 (the bath frequency), = 2iik/f3 (for k > 1) are Matsubara 
frequencies, c = ^(-i + cot £7/2), c k = 2^0^/(3(^1 ~ 7o)> fc > and = <4 + «m- 

In this equation, only po represents the physical system density operator, and other p n 
are called auxiliary density operators (ADOs). The subscript n is a multi-index, which has 
2 x (K+l) dimensions and can be extended as n m k, i.e. n w , nii,...,niK] n 2 o, n 2 i ) ... 1 n2K- The 
notation n mk ± 1 refers to an increase and decrease of this index. When the auxiliary density 
operators are all zero, the system and bath are not coupled; when the density operators are 
nonzero, we should properly take the coupling into consideration, which will influence the 
entanglement evolution dramatically. All of these values contain important information 
about the coupling. 

The method can be employed to calculate the dynamics of a system in a bath from any 
initial state (correlated or not) in a computer program with the truncating method developed 



by A. Ishizaki and Y. Tanimura [23|]. In the present work, we investigate the behavior of 
the entanglement death and rebirth time points with respect to several influential factors 
mentioned above. However, the range of all parameters should be carefully determined so as 
to obey the termination approximation and ensure that the results are sufficiently precise. 
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III. THE RESULT AND DISCUSSION 



We choose e = 1.5£, r] from 0A( to 0.8£, 7 from 0.4£ to 1.1C, and (3( from 2 to 3. Our 
initial state is chosen as 

/o o\ 

0.5 -0.5 
-0.5 0.5 
\0 0/ 

in the standard product |eie 2 ), \e1g2), | f7ie 2 ), \g1g2), and e m and g m are the excited state 
and ground state of the mth qubit. 



l*o (0) 



(9) 



A. System-Bath Interaction 

We know that rj is related to the system-bath coupling strength. Though the function of 
f] in the entanglement dynamics we can infer that when the coupling is weak, a slow decay 
of the concurrence is expected because the flow of information from the system into the bath 
will be slow, whereas the decay will definitely be rapid when the system and bath are tightly 
correlated. Figured] proves our inference; we see that the decay of occurrence is slower when 
77 is small than when 11 is large. The figure also indicates that small 77 may result in strong 
vibration of the concurrence. The environment can return a small amount of information 
to the system during the evolution. With weak coupling (small 77), this little returned 



o 
o 



■T!=0.2? 
r,=0.8? 




FIG. 1: Concurrence as a function of time with /3£ = 2.5, 7 = 0.5£, r\ = 0.2£ for the black line and 
rj = 0.8£ for the red line. 
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FIG. 2: Entanglement death time and rebirth time as functions of rj, with 7 = 0.5£ for the the 
solid and dotted lines and 7 = 0.4£ for the dashed and dash-dotted lines, /3£ = 2.5. 

information will have a strong effect on the curve of the entanglement when the concurrence 
is small, and thus, there is oscillation. As rj increases, the curve becomes smoother, and 
there is only one pair of death and birth time points. From Fig. [5J we see clearly the change 
in death and birth time points with the increase in 7]. The negative slopes of these curves 
are consistent with the effect of rj discussed above. It is interesting to find that at some 
intervals the slopes of these curves suddenly change in Fig.[2j This phenomenon is due to a 
change in the slope near the death and birth points on the concurrence curve. 



B. Characteristic Frequency of the Bath 

FiguresE] and H] show the influence of the characteristic frequency of the bath on the 
dynamics of the entanglement. Large 7 means a fast bath. Because the frequency of the 
two level system in this report is 1.5£, and if we take into account the range of 7 discussed 
here we know that the bath is slow. Most theoretical treatments with which we describe 
this entanglement dynamics are valid only if there is a fast bath, which has a much larger 
characteristic frequency than the system. Only in the fast bath scenario can we choose the 
initial state in which the system and bath are not correlated IjJ. However, as mentioned 
before, the hierarchy method can easily solve a slow bath situation, taking into consideration 
the bath effect on the entanglement in the initial state. A slow bath can receive and return 
the information from the system slowly and a fast bath may absorb all the information in a 
very short time and there may not even be revival, as we show in Fig.[3j This can explain 
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the delay in the death and birth when 7 is small. More details are shown in Fig.Hl The 
death and birth time points approach zero with an increase in 7. With smaller rj, we see 
that all time points increase, for the reason discussed above. 




t(l/0 

FIG. 3: Concurrence as a function of time with (3( = 2.5, rj = 0.6C, 7 = 7£ for the dotted line and 
7 = 0.5£ for the solid line. 
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FIG. 4: Entanglement death time and rebirth time as functions of 7, with rj = 0.5£ for the dash 
and dash-dotted lines and 77 = 0.7C for the solid and dotted lines, (3C, = 2.5. 



C. Temperature 

Temperature influences the system dramatically. In Ref. [l^] it is seen that the low- 
temperature scenario is unsuitable for Born and ultrafast bath approximations, which are 
used in many theoretical treatments; however, this scenario is in the solvable range of the 
hierarchy equations approach. In Fig-EJ the death time point increases and birth time 

7 



point decreases with a decrease in temperature. That is to say, the entanglement revives 
more quickly at lower temperature. The figure infers that the concurrence curve should be 
sufficiently smooth, and with the increase in (3(, the whole curve is raised up to obtain a 
smaller death time interval. If we simulate the two curves in the figure and calculate the 
relationship between f3( and the time point, then we may simply change the environment 
temperature to prepare our desired states and dynamics of the entanglement in experiments. 

In fact, with higher hierarchy of the ADOs, we can investigate a larger f3( scenario. 
However, the computation time is too great as it increases exponentially with higher orders. 
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FIG. 5: Sudden death time point (solid line) and sudden birth time point (dotted line) of entan- 
glement with different values of (3( from 2 to 3. Other parameters are rj = 0.6(, 7 = 0.5£- 



IV. CONCLUSION 



In this report we studied the entanglement dynamics of two-qubit in a non-Markovian en- 
vironment using the recently developed hierarchy equations approach. We explored in detail 
the role of parameters in the entanglement evolution, including the system-bath interaction 
77, the temperature (3 and the bath frequency 7. All these important features are useful when 
preparing a desired system state in an experiment. We also showed this hierarchy equations 
approach to be effective in calculating the dynamics of a system in many possible scenarios. 
Applying a higher order of the equation for a broader range of parameters and the physical 
meaning of non-zero ADOs will be future work. 
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